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Introduction 




This note describes two computer programs which are available for 



A 

h 



factor analysis using the General Electric Time-Sharing- Computer System. 

One program obtains the principal components of a correlation matrix; the 
ither makes an orthogonal transformation of the principal components ^ . 

using a varimax algorithm. The primary purpose of this paper is to facxlitate 
use of the programs ; for theoretical interpretations and mathematical 
detail the reader should refer to the books listed in the bibliography 
section. With regard to computational aspects, some familxarxty wxth the 



Time-Sharing Computer System is assumed. 

.. The factor analysis package consists of two completely separate 

programs named PCA (Principal Components ^alysis) and VARR (Varimax Ro- 
tation). In most applications the analysis would proceed according to the 

following sequence: 1. Run PCA. 2. Vi^SliX. 

and choose the principal components to be rotated, '3. Run VARR. Either 

program, however, may be used alone if desxred. 

The basic input to PCA is the correlation matrix which summaries 

the original data. The terminal output includes the principal component 
vectors, listed in order of the variance accounted for, and the latent 
roots of the correlation matrix. Some additional output is optionally 
available. PCA also' writes the principal component vectors (referred to 
collectively as the principal component matrix) in a permanent file named 
-PC (Principal Components) where they may be used later by the program 




VARR. 






The basic input to VARR is a set. of principal component vectors. 

The terminal output includes the rotated factor matrix, the percent of 

* « 

variance which is accounted for by each vector and some additional optional 

output. Jt 

V 

Prill cxpsl Comp 011611 ts Analysis 

There are a number of possible procedures for obtaining the 
principal components of a correlation matrix. The particular algorithm ■ 
followed by PCA is one which Horsti^ calls the Basic Structure Successive 
Factor Method With Residual Matrices. The method is a step-wise procedure 
which solves for one principal component factor at a time. The effect 
of, a given factor' is removed from the correlation matrix leaving a residual 
matrix from which the next factor is then obtained. PCA continues to 
eitract factors from the residual matrices until a latent root less than 
one is obtained; the factor associated with that root then becomes the 






last principal component. 

The procedure is iterative with respect to the solution for each 
principal component; That is, the solution begins with an initial 
approximation vector and proceeds through successive iterations with 
each cycle yielding a better approximation to the correct factor. Because 



/ P. Horst, Factor Analysis of Data Matrices. , (Neij York: Holt, Rinehart 

and Winston, 1965, pp. 156-167). 
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the procedure is Iterative, information must be supplied to the program 
to control termination. The detailed requirements are discussed in the 

gection on input and output . • — 

A 

'• An exact basic structure solution yields factors which are ^ 

mutually orthogonal to one another. The iterative method of PCA gives 
factors which are approximately orthogonal if enough iterations are 

carried out. 

The first principal component obtained by PCA tends to account 
for the maximum variance which can be accounted for by a single factor. 

If enough iterations are carried out, each successive factor accounts 

for the maximum variance not previously accounted for. 

■ ■ Other properties of the solution and an outline of the procedure 

are given by Horst-^. The purely computational parts of the program are 
almost identical to those proVided by Horst^^ The dimension and input/ 
output statements have, of course, been changed to conform to Time-Sharing 
FORTRAN. A listing of the program may be found in Appendix A. ^ 

Varimax Rotation 

Varimax is .the name given to one analytical procedure for making 
an orthogonal transformation. of principal component factors. As. a con- 
sequence of the varimax rotation the number of large and small loadings 



2/ Ibid., pp. 157-167. 

3 / Ibid., Appendix, p.607. 
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tend to be maximized and the number pf intermediate loadings tend to be 

■ minimised. In the foregoing sense the final factor matrix is regarded 

* * ' ’ 

as a simpls structura. 

There are several possible variations of the varimax computational 

4/ , 

=' ' procedure. ' The one used in VAER is referred to by llorstr- as the 

Successive Factor Varimax Solution. The new factor vectors are obtained 
in a step-wise fashion. The first factor tends to be that one for which 
the variance of the squared loadings is maximized. Subsequent factors 
are mutually orthogonal and, subject to this constraint, tend to come 
out in the order of the variance of their squared loadings. 

' The procedure followed in VAER is the so-called normal varimax 
in that each row of the principal components matrix is divided by the 
corresponding commonality. The final varimax rotated matrix is 
thus normalized. Ifshould be noted, hcxvever, that the Successive Factor 
Varimax Solution does not give the same answer as the varimax method 

originally suggested by Kaiser-^^ 

As in PCA, each vector is obtained by an iterative procedure . 

and information must be' supplied to the program to control truncation. 

Details will be found in the input/output section. 



S'FfKkre;,'’'^omputer Program for Varimax Rotation 
“ Educational and Psycholog ical Measuremejit, Vol. .19, ho. 3, Aut 

pp, 413-420. 
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Restrictions 






The principal components analysis will accept up to fifty variables 
in the correlation matrix and will coir>ute up to thirty prihcipal components 
/'Extraction of principal components is automatically stopped, however. ^ 
after appearance of the first latent root less than one. 

the varimax analysis will accept up to ten principal components 



.n 



with up to fifty variables. 

Input and Output 

PCA Program 

Part of the input to PCA is via a data file named EMAT (R Matrix) 

and part is entered on-line during execution of the program. 

• ■ Prior to running PCA, the upper triangle of a correlation matrix 

is entered in EMAT in Time-Sharing FORTRAN Standard input format. The 
contents of a sample EMAT file are shown in Appendix B. Only the sequence 
of values is important; the arrangement of data with respect to line numbers 



Is optional. ” • 

After the correlation matrix has been entered in RMAT, PCA is 

ready to RUN. During execution of PCA a request for five more items of 

information will appear on the teletypewriter in the following fprm: 
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#■ C*'*' ''Tri>^'t, 
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N > I NMAT> I f^£S j PjLI- 

•k . %• ♦ 

• * ■ ^ 

# 

# 

{ The user must then enter five numerical values from the key-_ 

r 

f\' . 

^ board. The symbols have the following meanings. ^ 

N The order of the correlation matrix entered 

in KMAT. 

INMAT A programming variable which controls the output. 

If INMAT=1, the correlation matrix is printed 
as part of the output; if INMAT=0, printing of 
the correlation matrix is suppressed. 

IRES A programming variable which controls the output. 

If IRES=1 (and INMAT=1) , each residual matrix 
Is printed out; if IRES=0, the residual matrices 
are not printed. If IRES has a value of one then 
so must INMAT. 

■ p ■ .A stabilization limit on the number. of iterations * 

perform.ed to extract any given, principal component. 
After each iterative cycle a "measure of improve- 
ment**^^ is computed and compared to the value of P. 



.'6/ For a details regarding the "measure of improvement" see the 
*”• listing 9 f PCA in Appendix A or Horst, pp«ci^ 



When the"measure" is less than the value of P , 
the iterative procedure^ is truncated. The user • 
may enter a positive value if he chooses or by 
; entering a zero the program automatically sets P 

equal to .00001, a value which will probably be 
satisfactory for most purposes. 

LI An absolute limit on the number of iterations 

performed to extract any given principal component. 

The user may enter a positive integer if he 
chooses or by entering a zero the program 
vill automatically set LI equal to 30. 

After entering the required data the printout at the teletypewriter 
might then appear as follows (user supplied information is underlined) : 

\ N » rMMAT# I P.» LI = 




The output of PCA is as follows : 

CORRELATION MATRIX. (optional). The lower triangle of the input correlation 
matrix. 

PRINCIPAL COMPONENT i. The i^^ principal component written as a row 
•vector. The loading on the first variable appears first and so on. 

i 
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RESIDUAL MATRIX i. (Optional) The residual matrix after the i prin- 
cipal component has been extracted. .Written in lower triangular form.- 
NUMBER OF ITERATIONS. The number of iterations required to extract 

i each principal component. Written in the order of extraction. 

\ 

i LATENT ROOTS. The latent root associated with each principal v- 

component. Written in order of extraction. 

CUMULATIVE PERCENT OF VARIANCE ACCOUNTED FOR BY COMPONENTS. The 
cumulative percent of total data variance which is accounted by the 
first i components. 

See Appendix C for a sample output. 

• *■.. PC Data File 

f 

The principal components extracted by PCA are stored in a 
file called PC for later use as input to VARR. The components are 
stored in Time-Sharing FORTRAN standard format in the order in which 
they were extracted, i.e.‘ , the first component extracted is first 
. in the file and so on. Contents of a sample PC file may be found' in 

Appendix C. . « 

Prior to running PCA, a dummy PC file must be set up large 
enough to contain the principal components which are output from PCA. 
The Time-Sharing Library contains items of six different sizes which 
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may be used to set up the dummy file/ 

A set of principal components. may also be entered directly 

Into PC rather than as output from PCA. 



VARR Program 

Part of the input (viz. , the principal components) to VARR 
is from the file PC and part is entered on-line during executxon of 



A 



the program. 

A set of principal components must be available in PC prior 
to running VARR. Procedures pertaining to PC have already been 



discussed. 

During execution of VAEK' a request for six more items of information 
willappear on the teletypafriter in the following form: 



ITKANS> P>ML = 



o . - 



7 General Electric Company, Time-Shari ng FO RTRAN Reference Manu^ ,, 
(General Electric Company, October 1966, pp. 65-66^. 
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The user must then enter six numerical values from the key- 
board. The symbols have the following meanings: 

M The number of principal components to be 

rotated. The first M principal components 
.In PC will be rotated. 

N The number of variables represented in each 

principal component. 

INMA.T A programming variable which controls the 

output. If INMAT=1, the M principal components 
are printed out; if imT=0', printing is 



suppressed. ?•. 

* ITRANS A programming variable which controls 

the output. If ITRANS=1, the final trans- 
formation matrix is printed out; if ITRANS=0, 
printing is suppressed. 

. p A stabilization limit on the number of iterations ^ 

performed to obtain any given rotated vector. 

After each iterative cycle a "measure of 
Improvement" is computed and compared to 
the value of P. When the "measure" is less 
than the value of P, the iterative procedure 



For details regarding the "measure of improvement" see the listing of 
VARR in Appendix D or Horst, op .cit ^ 



isttuncated. The user may enter a p 9 sitive 
value- if he chooses or by entering a zero 
the program automatically sets P equal to .0001, a 
value which will probably be satisfactory for 

% 

A 

most purposes. 

An absolute limit on the number of iterations 
performed to rotate any given principal 
component. The user may enter a positive integer 
if he chooses or by entering a zero the program 



will automatically set NL equal to 40. 

• After entering the required data the printout at the teletypewriter 
might then appear as follows (user ^supplied information is underlined) : 



K , I I T RAMS> ?* NL= 



The output of VARR is as follows : 

principal component matrix, (optional) . The M principal components 
selected for rotation. Each component is in column vector form. 
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NUMBER OF ITERATIONS. The number of iterations required to obtain 

each rotated factor. Written in order of their transformation. 

« • 

PERCENT VARIANCE ACCOUNTED FOR BY EACH FACTOR. The percent variance 
accounted for by each factor with respect to only the variance 
accounted ior by the M principal components. That is, the sum of 
these percentages add to 100% neglecting round-off error. 

FINAL TRANSFORMATION MATRIX. (Optional) The matrix which transforms 
the principal components matrix to the varimax rotated factor matrix 
Written as a series of column vectors. 

VARIMAX ROTATED FACTOR MATRIX. The new factors written as column 
vectors in the order of their transformation. 
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APPENDIX A 
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T.1 stin g of PGA Program 

100 DIfiE.'iSIOM rU 1275)^ A(5.-3)^ 

110 SFILS i^xMAT^PC 
1 20 PRI MT'*i\^.I I R‘£S^ L I =” 

1 30 I M?IJT.» I -X'l'.AT^ I R£3^ P^ L I 

• 1 40 IFC ?) 6 /■ 

/ 1 50 5 P=. 60001 
-,1 60. 6 IFCLI) 7^7^8 
/ ’ 1 70* 7 LI =36 ■ ' 

180 8 . M=(N'j=CM+l> )/2 . - 

^ 190 RFAD Cl) ( R( I ) ^ I = 1 ^ M ) . 

2 06 £=0. 

2 li^ DO 562 L=1^30 ‘ 

2 20 ■ IFCIi'^MAT) .10^10 • ' 

236 1FCL~1)‘ 86^80^9 • 

• 2 40 86 PRINTt."COR21£LATION MATRIX" 

• 250 GO TO 16 

260 9 IFC IRES) 10:i 16 

276 PRIMTt"R£SIDlIAL MATRIX"^ 

2 80 PiRIiMTj L~ 1 

■ 2:90 10 LK=L' -V . 

'3 0 0 D 0 1 0 2 I — 1 ^ M . ■ ^ . , 

3 16 DO lul J=l.» I ■ . 

320 IJ=ItC ( J“ 1 )=?=(^^•'=2•rJ) ) /2 
330 101 UCJ) = RCU)' 

3 40 . IF ( I X’MAT ) 1 62 ^ 1 02 ^ . ■ 

350 IFCL-l)' 1020>1020‘ 

360 IFC IRES) 162^102 ' 

■ 370. 1020 PRINTS tCUCJ)^J.= l^I) 

: 380 102 CONTINUE ' • ■ . 

390 DO 18 I = 1/M ' ' . ; 

■ 4 00 13 ACI) = 1 . . ■ ' ■ ' . 

4 10 DO 41 K=1^LI ' . ■ 

420. Op 28 I = 1^N 
430 11=1+1 

440 UCI>=0. . . ' 

• 4 50 DO 24 J=l^ I 

460 IJ=1 + C C J-l)’i=CM=5^2-J) )/2 
470 24 LICI)=JCI) + RCXJ)*ACJ) 

4 80 DO 28 J=IUN 

490 IJ=CCI“l)>j^C M V2 - i ) ) /2+ J 
560 23 iJCI)= UC I)+RCIJ)^A.C J) 

5 10 S=0. 

•520 DO 31 1 = 1 
-530j31 S=3+'JC I ) i*-AC I ) 

5 40 S= 1 • /S QRTF C S ) 

550 DO 34 I=1^X , 

566 34 ACI)=J(i)=-:^S . ‘ 

5 70- 3=0* 

586 DO 37 1=1 

5 9 u 37 3= S + A C I ) ^- ❖ 2 ; . 






o 

ERIC 



I 






I ! 

&. 



^ ■ 

: 

i - 



6 00 

6 10, 
629 



^630 

i640 



/. 



A •6 50 



jt 



. 6 60 
L 6 79’ 
6 80 - 
690 
7 99 

• 7 10 
7 29 
7 30 
7 40 
■ 7-50 

■ 7 60 
770 
•. 7 80 

7 90 

8 00 
'8 19 

820 
8 30 

• 8 40 
8 50 
8 60 
870 
880 
890 






E.1=F. 

E=SQRTFCS) ^ /I 

IFC (A83FC£l/ii:**l 

41 COMtIi''jJ£ . 

kvcl)=k-i ■ : 

GO TO 43 ^ 

42 K\/CL)=K ' 

43 EVCL.)=S 

P.RINTf'PRIMClPAL COWPOMEM . 

PRINTS L , • 

PRINTS t C A'C I ) .» I = 1 ^ ^ 

WRITE C2) CAC I >,» l = l-»i'^^ 

IF CE-T"*5 52jA7j;A7- 
47 00 50 1 = 1 
DO 50 J=1-»I 

IJ = 1 4 - ( ( J “ U C iM*2- J ) ) /2 

50- rciJ) = R<iJ^-:^^^^"'^^'J^ 

IFC 1 • ~S) 502,» 52,»’52 
50^ CO(vYIi^U£ 

52 pRIiN5Tt"i'i'JF'P£P ITERATIONS' 

PRINT,» t (KVC I ) •» 1 = 1 •. • • 

PRI NT t "LATENT ROOTS" . ■ v. . . 

PRINt,» f CEVC !),» I = 1->TK) ••.• 

PERC 1) = 1 90 • ^ ^ ■ ■ ■ ■ - ’ 

bo 60 I=2^LK . . 

PRINT,» ^ CPERC !),» I-1-»LK) ' • . •. 

STOP .. ' ■ ■ . '• . ■• .■ 

ENo . ' . , • • • • •. • ■ ■ ■ ■■■ • ; « 
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APPENDIX 'B 






Sample PCA Problem 
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The sample problem is one given by Horstr" . The correlation 

matrix is : 

1.000 






Ji 



.829 


.768 


.108 


.033 


.108 


.298 


.309 


.351 


1.000 


.775 


.115 


.061 


.125 


.323 


.347 


.369 




1.000 


.272 


.205 


.238 


.296 


.271 


.385 






1.000 


.636 


.626 


.249 


.183 


.369 








1.000 


.709 


.138 


.091 


.254 








% 


1.000 


.190 


.103 


.291 












1.000 


.654 


.527 



1.000 



.541 

1.000 



A listing of the input as it appears in the data file KMAT is 
as follows : 

_• « 

■ 1 10 1 . .82-9^ .768^ .163^. 03*3^ . 103^ .298^ .309^ .351^ 1 «775 

.,120 .1 15^.061^. 125^ .323j.3A7j .369^ 1 .000,» .272^.205^ .233^ 
i 30 •21\* • 385.» 1 • 06O.» . 636j • 626j »2A9j • 133.» • 369> 1 • OG-0.» .709 

1 40 • 133.» .091/. 254> 1 • 000.? • 1.90.? • 1 03.? . 29 1 ? 1 • G00.» • 654.? •527 j 1 • 

1 50 *.541> 1 ' ' 



■ 9J Horst, Ibid., p.l22. 
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The values for N, INMAT, IRES, P and LI are 9,1, 1,0,0, 



respectively and are supplied on-line, 



A 



1 



1 



;’5 

‘•1 






n. 



■ i 



, 



,.j 



■^1 



■ % 



■ I 






r .JV' 



■I 
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APPENDIX C 



Sample PCA Output 

Printed below is the teletypewriter output from a principal 
components analysis of the problem given in Appendix B. Information 
■supplied on-line by the user is underlined. 



A 

V 



. PCA 14:03 



N ^ i AT ;»IR£S>P;»LI = 
?• 9.> 

CORf^HILATIOW MATRIX 
' . ■ 1*00 

. v829. . 

^ • * 

-.768 

• ^.108 

• ' ■ . .033 

-.108 

1.00 

.298 . . 

.19 

.309. 

.103 

.351 

.291 



W1 THU 10/05/67 



00 


t 


775 


1.00 

r 


• • 

115 


.272 


\ 

06 1- 


■ -.205 


125 


.238 


323 


..296- 


00 




347 


.271 


654 


1.00 


369 


' \ 

.335 


527 


^541 



1 i 00 


t 


.636 ■ 


• ^ 1.00 


.626 


■ .709 


.249 


.133 


• 


1 


\ 




'. 183 


.091 

• 


■' .369 . 


.• 2d4 


1 . 80 


• • > 



P RING I PAL C OMPONEiN'T 

, f 


1 


: . .7152 


.7335 


.5196 


.6403 


' * \ 


.* 


RESIDUAL MATRIX 


1 ■ ' 


.4385’ 


• 


53008 ■ 


.4546 


52159 ■ 


.2049 


- . 29-06 ■ 


-.2966 

! 


■ .-.2993 


-.2827 


-.2636 


-.2538 


,.73 




/ ■ 


• 


■' - --.1599 


-.1499 


; , -. 1427 


.59 


-.1305 


-.1068 


-.2163 


.2605 


-.1604 


-.1591 


-..0306 


.0691 


■'P RI MCIPAL COMPONENT. 


2 


■ -.4965 , 


-.4318 ■ 


. . 69 12 . * 


-.0784 • 


. RESIDUAL MATRIX 


2 . 


• '.*242 




.0616" 


.2225 


• ' .067 


V. *' 

-.0604 


• • ■ .0306 


.0151 


'1-0685 ■ 


. .0747 


.0795 . 


. .0743 


.2522 




-.1939 


-.1877 


■-.0o35 


.5833 


-.2121 


’ -.1861 


-.1026. 


.2476. 


-.1433 


-.1,43 


,-.1033 


.0717 



.772 .557/^ 

.6145 .7151- 



..4041 • ■ ' 

-.15S3 .6893 

A 

itr 

-.15*42 .3766 

.-.1631 ■ .33^64 . 



T. 1983 ■ 


; -. 1079 


-.2034 . 


-.159.5 


. 6224 


. ^ 


-.1671 


-.0296 


.1016 


. .4386 


-.30 


.647 


-.1645 

/ , 


. .0335 




* • 






.3141 




,.0358 


-.2708 


.8683 


--.1033 




■ . 1 ,: 


.0442 


-.1103 


-.2213 


-.0572 


-.2527 


-.053 


.5954 




-.1.57 


-.0513 


.1071 


.4375 



.4654 



.7834 

# 

. 4672 
.16 
'.1949 
.0788 

I 

.7413 



.2332 

.0455 

. 1013 

.0729 

. 1-037 



V 



P 21 -NCI PAL COMPONENT 


3 


• 




-.3478 

: -.1918, .. . 

• 


-.319 
. 5883 


-.4643 . 
.6222 


-.071 • 
.3677 


RESibU.AL MATRIX 


3 * 


f 




‘*.1211 


• 


f 

• . . ' 




-.0493 


. 1207 
/ 






' -.0733 


/ 

-.0688 


. 1502 

* . • • • • • 




1 

. • .0.059 ’ . 


-.0076 


. 007 ' 


.2657 


• . . 0043 ■ 


'.0159 


. -.0064 


-.1164 


.0123 


.0131 


-.0335 •. 


-.■ 1244 


'.2155 






. 


.6057 2.61904E-05 


.0163 


-.0154 


.0243 . 


.2378 


■ / ^ 


• 


• \ 

^ • .0042 

• .0168 


.0124 - 

-.1 184 

■ / 


8. 18963E-04 

-.2082 

4 


-.0089 
^ • 


t » \ , 

“ . .... ■ -.0159 

• \ -.0333 


-.0257 

-.1445 


■ -'.0082 
-.1216 


• -.0252 
.3523 


PRINCIPAL COMPONENT 


\ 4 


• 


t 


. .0305 
. . . 1312 


.057 , 
.2’845 


-.0094 

.1661 


-. 1074 
-.5648 


iNI UMBER OF ITERATIONS 


: - * - * 


; 


• " * 

7 


8 


8 
■ "• 


16 


.latent roots 




• 




, '3.7491 


2.0495 


1.3308 


■.’47 44 




A 



€ 




CUMULATIVE PEiWENT VARIANCE ACCOUNTED FOR BY COMPONENTS 
4i.6562- 64.4286 79.2154 84.4866 







•AT’ LIi\’E NO* -86G: STOP 



:,RAN 172/6 S£C. 



; 1844 



. ./ • ■ 
.1992 

.0809 
.0067 
.0419 ' 
.0359 

. 1175 
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A 


listing of the 


data 


file PC 


as written by the 


program PCA 




appears 


below. 






• 






1 000 


.7152 




.7385 


.772 


.5574 


.4654 


.1010 


.5196 




.6403 


.6145- 


.7151 




1020 
•% • 


-.4965 




-.4818 


- . 3.0 


. ...647 


..J418 


1030 


.6912 


• 


-.0784 


-.1645. 


.0335 


• . 


/l 0^0 


-'.347S 




-.319 


-.4048 


-.071 


-.1844 


1 056 


-.1918 


■ 


. 588 3 


.6222 • 


.3677 


* * 


1G60 


.0305 




.057 


. *:.G094 


-. 1074 


.1175 


1070 


.1312 




.2845 


1661 


-.56.48 








t. 







■i 



m 



> 




APPENDIX D 



Listing of VARR Program 

# 

Several comment lines in the listing, below have only numbers 



on them. These numbers may be used to relate the subsequent FORTRAN * 
statements to the numbered equations in Horst ^s text— For 
example, line 460 has the number 1; this means that the next few 



statements pertain to equation 18.3.1 of the text. 



1 00 DIMENSION -AC50^ 10)^ BC50^ 10) 
110 +aS( 10) ^ VK 1 0) ^ UC.l 0) > PERC 10) 
1 20 SPILE PC 



,x4(! 10^ i0)^ C(50)^ V(50)^KV( 10) 



i 30 PRI I I TRANS^ Pj NL=" 

i 40 I NPUT^ iAj I NMAT^ I TRAi\'S^ Pj ^’L 

1 50 IFCP) 10^ 10> 11’ 



160 10 P=.0001 



170 11 IPCiX'L) 12>12^13 



1 80 12 NL=40 

190 13 READCl) CCAC J)^ I = 1^M)^ J=1^.M) 



2 00- IFCINi'iAT) 25^25 - ' 

2 10 PRIi\jYt’'PRIiNiCiPAL COMPOME^'T MATPvIX” 

-220 JA=1 ■ • 

230 IFCM-3) 20^20 . .• 

2-40 J3 = 5 • ■ ^ 

250 GO TO ’21 ' ^ , - « • 

.2 60 20 JB=M . . . • - . 

270 21 DO 22 I = 1^M ’ ’ 

280 ’22 CACi’^ J)^ J=JA^ JB) ’ ' 



290 IFCJB-Fi) 23^25^25 
.300 23 PRINTS t 
310 JA=J3+1 - 

320 JB=JB-f-5 
3 30 IFCJB-M)- 21 
340 JB=M 
350 GO TO 21 



3 60 * ROl-/ MOR^IALI BATI ON 
370 .25 00 1 15 1 = 1.»N . ' 



10/ P. Horst, Ibid, pp. 424-425. 






' 3.S0’ C( 1 

■ 390 DO 1 15 

40‘0 115 C(1)=C(I)+A(I^ 

■ 410 DO 117. I = 1^N 

420 117 'C(I) = l«/30ivTF(CCi)) 

430 DO 120 I=1^N 

440 DO 120 J=UM ’ ' 

4 50 120 A(I^ J^=0^ J) . 

460* 1 . ‘ ' •• 

.470 DO 127 0=1^^!’ * ’• 

-480 W(J)=0* 

;'490 DO 127 1 = 1^ N 
5’00 1 27 vj ( J ) =VJC J ) +A( I ^ ,J ) 



5 10 FiN=N . . 

' 520 AL=0. 

• -530 DO 131 J=l^ M 

540 131 AL=AL+’'K J) 'j; ‘i^2 
550 AL=1./SQRTF(AL) 

560 DO 134 J=1^^5 
570‘ 134 WC J) = '.'K J)*AL 
•580*2. 

590 DO 138 I = l^-N 
600 VCI)=0« 

610 DO 133 J=1^M 
- 620 i38 \/<I)=VCI)•^•A(I^ 





•630*3 . . ' ' 

640 DO 1-78 L=l^Fi’ • . . 

•650 LA=1 

6 60 DO 144 1=2^ iM \ 

670 IFCV(LA)-7(I>> 144^144^143 

680 143 LA=I 
,•690 144 COMTIMUE 
■ • 7 00*4 • • 

7 10 DO 146 J=1^M 

720 146 H( L)=A(LA^ J> 

7 30 £=0* 

; 7 40 DO 173 K=l,»i\’L 
7 50 KV(’L)=K ' 

7 60*5 . 

770 DO 153 I=1^N 
7 80 ,B( L.»L) = 0« 



790 DO 153 J=1^M 
8^;0 153 BCI,»L)=3 

810*6 ■ I 

820 8SCL)=0* 



C I^D+ACI^ J^L) 



830 DO 156 I=1^N 
'840 156 SS C L) -B3 ( L) ( 
.853 -BSCD^BSCD/FM 
860 £1=3 



I,L)' 4 =:i =2 



870 £=AL ^ • 

380 IF C P~A3SF ( £ 1 1 • ^ ^ 
8 90 158 DO 159 I = l^i''' 

•900 159' B( I.»L)=3C 



158^ 15S,» 174 
3 - 3 a^L)^i=BS(L) 



I 



1 



I' 



I- 



I 



k 






.} 

I?'- 



P: 

I" 



ti 






I 



i' 



I 






«: 






f' 



I 



er|c 
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/ 



. J 



9 16*7 

9 20 -DO 163 J=l^>i 
9 30 UCJ)=£. 

9 46 DO 163 1 = 1^ i\’ 

950. 163 U(J)=UC J)+8C I:.L)^AC 

960*6 , ' • • 

970 AL=a. . 

989 DO 166 J=l^:^ ‘ ■ ’ • 

9 90 166 AL=AL-HJ( J)'<=^=2 • • 

1000 AL=SQRTFCAL) 

101'0*9 ’ ■ ■ 

1 020 D.O .169 J=1^M ' 

“1039 169 HC J^L)=IKJ)/AL 
1040 173 COMTINUS .. ^ . 

1 050* 16 ■ . 

1060 PRINT 

1070 174 DO 176 i = l^N 
1080 DO 176 J=1^M . • ■ ' 

1690 176 AC I ^ 3 ) “ A C I ^ J ) ~B C I ^ L ) '*'H C «J^ L ) 

1100*17 ' 

1110 DO 17S I = 1^N . . . . • 

1 120 176 VC I)='yCI )+BC I^L) 

1 130. PRINT T‘*NLK-i3£R OF ITERATIONS” 

1146 PRINTS CKVCL)^L=1.»N) : ' 

1 150 DO 179 L = 1^M ^ . 

1 166 179 ■ PERCL) = 169.G=>BSCL) . 

.1 170 PRINTS "PERCENT VARIANCE ACCOUNTED FOR BY EACH FACTOk 

1180 PRINTS CPERCL)^L=1j;N) 

1190 IFCITRAM3) 190.»196 






A 



\ 



1200 PRINT t "FINAL TRANSFORMATION MATRIX" 



- 1216 JA=1 
1220 IFCM-'5)^ 184^134 
'1230 jO = 5 ' . 

1246 GO TO 185 . ^ 

1250 184 J3=M 
. 1260 135 DO 186 I=1^M 
1270 186 PRINTS CHCI^J)^J=JA^JB) 

1230 IFC JB-M> ,187.» 190 190 

1'290 187 PRINTS t 
1 300 JA=J3+1 
1316 JB=JB+5 
1326 IFCJ3-M) 135 
• 1330 JB=N 
1 340 GO TO 135 

135G 190 CALL SINCHAC A^ 8^ ’•/.» V^ N.»M.»NL^ P) 
1360 STOP 

1 37 0 END ' .... '.... ' . • • 







25 



13S<3 
1 390 
1 400- 
1410 
1 420 
•i 430 
1 440 
1 450 
•1 460 
1'470 
.1 480 
"a *490 
‘1 500 
1510 
1520 
1530 
1540 
1550 
1560 
, 1570 
1580 
1590 
1 600 
1610 
1 620 
. 1 630 
1 640 
1650 
. 1660 
.1 670 
1 680 
i 690 
i 



)) 12 ^ 12^15 



■1710 
1720 
1 730 
1740 
1 750 
1.760 
■1770 
1780 
1790 
1800 
1810 
1 820 
1 830 
1840 



SUBROUTI i\'E S I MCKA '( 

0Ii'*ili:M3I.0N A(50.» 10)^3(50^ 10)^ WC 1 
DO 6.I=1^M 
DO 4 

A(I^ J)=0* 

DO 2 K=1^M 

2 AC !.» J) =A( I J) 

4 ACJ^ I)=A(i:» J) . - 

6 VCI)=1. 

E=0. 

DO 14 L=l.»i\'L ■ ... 

DO 6 I=1^M • : , 

vJCI)=0. / 

DO S J=1^M '• 

8 ' IK I ) =vJC I ) +AC I .» J) J) 

AL=0.' 

DO 10- 1 = 1 

10 AL=AL+'KI)*V(I) • . 

AL=30RTF(AL) ' . 

El=£ 

£=AL 

IF(P-ABSF(£1/E-1 
12 DO 13 I=1>M 
13. VCI)=W(I)/AL 
14 CONTINUE 
15D016I=1>N 

u(i)=0. 

•DO. 16 J=l.»v''i 
16 VKI)=v/(I)-i*B(I.» 

DO is 1 = 1^ N .•> 

DO IS J=1^N 

3=SIGNFC 1 )=5^SIGNF( 1 ..»'•/( J) 

.PRINT t.”VARIMAX PvOTATED F.ACTOR 
JA= 1 

IFCM-5) 20.»20 / ■ 

JB=5 ■ . ■ ■ , ;^..- 

GO TO 21 

20 JB=N 

21 DO 22 1 = 1 .»N 

22 PRINT.» (3( I.» J).» J=JA.» J8) 

IF(J3-rO 23V25>25 , . 

23 PRINTS 7 

JA=JB+1 ’ ■ • 

J3=JB+5 

IFCje-M)- 21 

JB=M 



NL.» 

0 ).» 



P) 

VC 10) 



MATRIX 



tt 



1850 GO . TO 21 
I860 25 RETURN 
1870 END 
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APPENDIX E 



Sample VARR Problem 

The problem for VARR is to rotate the .foliating principal 



4 components: 




• 






• ^^4 


• 1 
/ * 

, Principal 


Principal 


Principal 


Principal 


Jt 




Component 1 


Component 2 


Component 3 


Component 4 






.7152 


-.4965 


-.3478 


.0305 




• 


.7385 

.7720 


-.4818 

-.3000 


-.3190 

-.4048 


.0570 

-.0094 






.5574 


.6470 


-.0710 


•"’-.1074 




i 

1 


.4654 


.7418 


-.1844 


.1175 




{ 

t 


.5196 


.6912 


■4 

-.1918 


.1312 






.6403 


-.0784 


.5883 . 


,2845 






.6145 


-.1645 


.6222 


.1661 






.7151 


.0335 


.3677 


-.5648 






A listing 


of this input 


as it appears 


in the data file PC is given 




below. 

.. .. ‘ 
1000 


.7152. 


.7365 . 


- ' , * % - 

.772 


.5574 


• 46 


.1010 


.5196 


• 6403 


.6145 


.7151 


• 1 A 

\ 


1020- 


.4965 


-.4818 


-.30 


.647 


1 030 


.6912 . 


-.0784 


- . 1 645 


.0335 


-. 18 


1 040’ • 


.3476 


-.319 


-.4043 


-.071 


1 050 


. 1916 


.5883 


.6222 


.3677 


.11 


•1 060 


. 0305- 


. 057 


-.0094 ’ 


-.1074 


1 070 


.1312 


.2845 


. 1661 


-.‘5643 






• ....»• ... 






. • 





J 
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The values for M 
respectively ^ and are 

i 

t 

/ • 



N, INMAT, ITRANS, P, NL are 4, 9, 1,1, 0,0, 
supplied on-line. 



jt 



m 
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•. appendix' F 



Sample VARR Output 

Printed below is the output from a varimax rotation of the 
principal components given in Appendix E. Information supplied 



bn-Hne by the .user .is _ under lined . 

14il2 W1 THU 10/05/67 




iM ^ I NMAT^ I TRAHS^ Pj NL = 



.PRINCIPAL CONP( 
.7152 
.7385 
• • • .772 

■ . .5574 

.4654 
.5196 
;6403 
^ .6145 

.7151 

\ 

NUMBER OF ITER 
5 



NENT MATRIX ■ 
-.4965 . 
-.48 18 
-.30 
.647 
i74lS 
.6912. . 
-.0784 
-.1 645 
.0335 

iTIONS' 

. 5- 



-.3478 
> -.319 
-.4048 
-.071 
-.1844 
-.1918 
■ .5883 

.•6222, 

. .3677 

/ ; 




• .0305. 
■, .057 
■ -.0094 
-.1074' 
.1175 
.1312 
.2845 
. 1 66 1 
-.5648 




PERCENT VARIANCE 
34.1572 



ACCOUNTED FOR BY EACH 
34*9125 23.831 

I 



FACTOR 

7.0993 



FINAL TRANSFORMATION MATPvIX • 
•• .5054 -.7039 

•:.8359 • --.5393 

-.2011 -.4612 

• - '.0741 • ’ *0315 



.4673 

-.0999 

.8462 

.2358 



-.1754 
- . 822 1 
-.1757 
' .9685 



VARIMAX' 



ROTATED FACTOR MATRIX 



.0199 
.0414 
•.2388 
.9597 
.9983 
’.9924 
.1752 
. 066$ 
.2783 



.9942 

.9335 

.9671 

.0842 

.018 

.6955 

.2512 

.2647 

.363 



AT LINE NO.. 1366: 'STOP^ 



RAN 



.1031 

• 1 455 

.0498 . 

.1279 

.0167 

.9474 

.9496 

.9603 

.5175 



184/6 SEC* 



-.0254 
.6681 
’.0724 
.2357 
— . w 5 3 4 
-.0607 
-.0671 
.0581 
.7503 



